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Abstract 


This report describes the analysis of 3275 Geos -3 arcs of altimeter data 
containing 624670 frame averages. This data was adjusted to remove orbit 
error and altimeter bias in a primary adjustment and four regional adjustments. 
The root mean square crossover discrepancy was about ±55 cm after the adjust- 
ment. The adjusted altimeter data, now considered to give geoid undulations, 
was used to predict values at 1° intersections from which an oceanic geoid map, 
with predicted accuracies, was prepared at a two meter contour interval. This 
geoid was compared to the GEM 9 geoid over very long profiles to examine the 
long wavelength error in the altimeter geoid. At a wavelength of 13010 km the 
root mean squares difference was 57 cm. The altimeter geoid was also com- 
pared to altimeter geoids fixed by precise orbits. We found a root mean square 
difference of about 1 m with a systematic difference that implied the equatorial 
radius of the earth was 6378137 meters. 

The adjusted altimeter data was also used- to determine a total of 29479 
1° X 1° anomalies ( and undulations ) , 27466 of which had an accuracy of 15 mgals 
or better. Their average accuracy was 8 mgals. In addition, 957 5“^ mean 
anomaly and undulation values were computed. Eepresentative anomaly differ- 
ences with terrestrial estimates were 12 mgals for the 1° x 1° values and 7 mgals 
for the 5° values. Additional computations of point anomaly values were made 
to compare with ship data in the area of the Ninety East Ridge. There we found 
the altimetry anomalies followed quite well variations of the anomalies at 100 km 
wavelengths, and clearly showed correlation with bathymetry. * 
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Introduction 


This report is an extension of the analysis with Geos-3 altimeter data de- 
scribed in Rapp (1977a). This extension was primarily oriented to processing 
an additional set of Geos-3 altimetry received during January and February 
1978, As will be described in subsequent sections this new data was merged 
with the previous data and used to obtain a near global oceanic geoid map and a 
large number of l°x 1° mean anomalies. In addition a number of computations 
involving point gravity anomalies was carried out for comparisons with actual 
ship measurements. 


The Data Set 


The data used in our investigations is the intensive mode data where we 
use frame averages. (A frame average represents a mean over 2.048 seconds 
or 3.277 seconds.) In taking the original data supplied by NASA ( through the 
Wallops Flight Center), an editing procedure is used that has six different edit 
criteria.' These criteria are designed to remove most of the bad or poor alti- 
metry data. In our first analysis (ibid, 1977a), we accepted 1976 arcs containing 
419,294 frame measurements. The location of this edited data is shown in 
Figure 1, 

During January and February 1978 a set of approximately 70 tapes were 
received containing additional Geos-3 altimeter data. The location of this un- 
edited data is shown in Figure 2. 

Portions of the data were arcs that were just the continuation of previous 
arc segments. A procedure was developed to merge common, continuous arcs, 
of the old and new data sets whenever possible. After editing and merging, the 
altimeter data set available for further processing is shown in Figure 3. This 
data set now contained 3275 Geos-3 arcs with about 624670 frame averages. The 
same data set is shown in Figure 4 which also shows the location and numbering 
of the 5° equal area anomalies used in this study. 


Bias Removal and the Adjustment Process 

The altimeter data that is available in the edited form is subject to errors 
caused by a bias in the altimeter and orbit errors. These errors must be re- 
moved (as much as possible) before additional processing is carried out. To do 
this a procedure was described in Rummel and Rapp (1977) and Rapp (1977a) that 
combined a fit of the altimeter implied geoid undulations to a satellite derived 
geoid with crossover constraints to obtain parameters of an error polynomial. 
After these parameters are found for each arc, adjusted geoid heights can be 
found. ( We note here that formally we are dealing with sea surface heights. If 
we neglect sea surface topography geoid heights would be found.) 
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Figure 1. Location of Edited Arcs in First Altimeter Analysis (Rapp, 1977a). 
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Figure 2. Location of Additional Unedited Altimeter Arcs. 
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Figure 3. Location of Edited Geos-3 Altimeter Data and Areas of Regional Adjustments 






For this work a procedure almost identical to that described in the earlier 
work (ibid, 1977a) was carried out. We first made an adjustment of a primary 
network of altimeter arcs that were chosen for their arc length and representative 
global distribution. After these arcs were adjusted, four regional adjustments 
were performed enforcing the primary arcs in the region. 

The mathematics of the adjustment process are described in detail in Rapp 
(1977a) as well as the procedures for determining the crossover locations. The 
reference field used for the reference geoid was the GEM 9 ( Lerch et al. 1977) 
set of potential coefficients taken to degree 20. For arcs whose length were less 
than 18® (305 seconds) only a single bias term was solved for. For the longer 
arcs a bias term and a trend term were found. 

The location of the regional adjustments are given in Table 1 and shown in 
Figure 3. 


Table 1. Regional Adjustment Locations 


Name 

Geographic Limits 



New Calibration 

70 to 10 

260 to 314 

East Atlantic 

70 to 10 

314 to 34 

Africa - India 

j 30 to 10 

20 to 120 

L 10 to 70 

218 to 120 

Pacific 

f 70 to-70 

120 to 218 

t 70 to 10 

218 to 260 


Information on the adjustment of the primary and regional networks is given in 
Table 2. 


Table 2. Adjustment Statistics Related to the Primary 
and Regional Adjustments. 


Region 

Number of 

Crossover Discrepancies 

Arcs 

Unknowns 

Obser- 

vations 

Cross- 

overs 

a priori 

a posteriori 

Primary 

700 

1383 

263077 

10149 

± 7.81m 

± 0.60m 

New Calibration 

478 

800 

81822 

35403 

6.94 

0.47 

East Atlantic 

273 

. 464 

47065 

9144 

8.56 

0.56 

Africa - India 

588 

1074 

139772 

15977 

7.42 

0.50 

Pacific 

408 

707 

95947 

12041 

10.38 

0.60 


i 






We shoidd note here that the adjustment process is carried out in two 
stages. A preliminary adjustment is made and adjusted crossover discrepancies 
are examined. Cross points having discrepancies greater than about 3. 5 meters 
are deleted. In some cases a whole arc segment may be deleted. The adjust- 
ment is then repeated a second time to obtain the data given in Table 2. 

The crossover discrepancies after the adjustment, as shown in Table 2, 
are of the same magnitude as found in the first adjustment (Rapp, 1977a, Table 
10 ). 

The adjusted undulations from this new adjustment (with additional data) 
were compared to the corresponding values from the earlier adjustment; The 
results of these comparisons are given in Table 3. 


Table 3. Comparison of Adjusted Point Geoid Undulation of New and Old 
(Rapp, 1977a) Adjustments. 


Area 

Number of 

Mean 

RMS 

Maximum RMS 

Arcs 

Difference 

Difference 

Difference * 

Primary 

337 

0.06m 

±0.33m 

±2. 04m 

Calibration 

154 

-0.09 m 

± 0 . 29 m 

± 1.78 m 


* Along any one arc. 


We see that on the whole the differences between the adjustments is on the 
order of 30 cm ( RMS ) although over a few ares the differences may reach the 
2 meter level. 


The Altimeter Geoid 


Disregarding sea surface topography effects, we may regard the adjusted 
altimeter data to give us geoid undulations with respect to an ellipsoid of defined 
flattening ( 1/298.256 ) but whose equatorial radius is specifically undefined. It's 
conceptual definition is, however, the equatorial radius of the ellipsoid for which 
the global mean geoid undulation is zero. 

A global oceanic geoid ( or mean sea surface ) has been computed from the 
adjusted data using primarily the procedure described by Kearsley (1977). This 
procedure predicts a geoid undulation (and its accuracy) at grid intersections 
from the surrounding point altimeter undulations. The predictions were made 
using least squares prediction techniques using covariances from subroutine 
COVA (Tscheming and Rapp, 1974) and including the data noise (with a ±0.4 m 
contribution from the effect of errors in the GEM 9 potential coefficients on the 
adjusted geoid) in the process. The grid interval chosen was l^'x 1® using the 
5 closest altimeter points. Choosing the grid interval this large can result in 
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the loss of detailed features below the scale of 100 km. However, the cost of 
computing a global oceanic geoid at a finer scale is prohibitive for us. 

The grid values were contoured using a new contouring program developed 
by Siinkel (1979) that uses spline functions to obtain smoother contours than seen 
in Kearsley (1977). The computations were made, roughly, in 30° x 30° blocks. 
All grid undulations and their accuracies were printed out and are available for 
other use. The result of the contouring (for the undulations and the accuracies) 
were output from a small ( 11" wide ) Versatec plotter. The individual blocks 
were fitted together to produce four map sheets (undulation and accuracy maps 
for the Eastern and Western Hemispheres ) whose original size was approximately 
39" X 28". Reduced versions of these maps are shown in Figures 5, 6, 7, and 8. 
Almost full size prints of these maps have also been made. Areas that are blank 
on these maps did not have sufficient altimeter data for preparing a geoid plot. 

No attempt here is made to interpret this geoid. The features obvious 
from potential coefficient geoids is present as are the details from short wave- 
length geophysical structures. 


Geoid Comparisons 


It is of interest to compare the geoid given in Figures 5 and 6 ( or more 
precisely the predicted grid values) with other sources. In the past comparisons 
with the altimeter geoid and some other reference geoid have been made over the 
relatively short arc segments of Geos-3 data. Our first interest here is to com- 
pare the GEM 9 undulations ( to degree 20 ) to the altimeter derived undulations 
in long profiles having a constant latitude or a constant longitude. For a number 
of different tests 7 long (-^14000 km) profiles were chosen. A plot of the GEM 9 
geoid and the adjusted altimeter geoid for two typical profiles ( two and five) 
are shown in Figures 9 and 10. ( The origin for profile two is at X = 10° , while 
the origin for profile five is at tp = -65°. ) Information on these profiles is given 
in Table 4. 


Table 4. Comparison of GEM 9 and Altimeter Geoid Profiles. 


Profile: 

Latitude 

Longitude 

Length 

Mean Diff. 

RMS Diff. 

Two 

Five 

-48° 

-65° to 52° 

10° to 186° 
335° 

14732 km 
14455 km 

-0. 6 m 
0.7m 

±2. 3m 
±2. 6m 


These plots do not reveal any significant systematic differences between the two 
undulation sets. 

The undulation differences ( GEM 9 minus altimeter derived) were analyzed 
by Eren (1979, private communication) to determine their power spectrum. The 
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Figure 9. The GEM 9 and Altimeter Undulation Along 
the -48® Parallel (Profile Two) . 
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Figure 10. The GEM 9 and Altimeter Undulation Along 
the 335° Meridian (Profile Three). 
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power, by wavelength is given in Table 5 based on averaging the results for the 
7 individual profiles. 


Table 5. GEM 9 Minus Altimeter Geoids by Wavelength. 


Wavelength 

(km) 

Approx. S.H. 
Degree 

Power 

(m^) 


GEM 9 
Error (m) 

13010 

3 

.33 

±,57 

O 

-H 

6505 

6 

.61 

.78 

.16 

4337 

9 

.60 

.77 

.38 

3253 

12 

. 58 

.76 

.43 

2602 

15 

.49 

.70 

.51 

2168 

19 

.15 

.39 

.53 

Total to S.H. Degree 19 

2. 76 

1.7 

1.7 m* 


* Due to errors in the GEM 9 potential coefficients from n = 2 to n = 19. 


The above computation assumed that the data was regarded as periodic. An 
alternate procedure assuming non-periodic data, with a window function was 
also carried out with similar results. 

From Table 5 we see that at long wavelengths we can expect errors in our 
altimeter geoid up to about 0. 75 meters. At shorter wavelengths (2168km) the 
difference is similar to the expected error. 

Additional computations could be carried out in this- area by using more 
profiles and introducing data noise. 

Our altimeter geoid was also compared to altimeter geoids presented by 
Brace (1977) and by Marsh et al (1978). The Marsh et al geoid (or mean sea 
surface) is based on the computation of certain laser reference orbits which are 
used to obtain reference altimeter undulations. These latter undulations are then 
used as a frame for a crossing arc adjustment process. The comparisons were 
made for Ihe three areas in Marsh et al where the plotted contour interval was 
1 meter. The values at grid intersections were read from the contour maps of 
Marsh, and taken as the predicted values from our data. The results are 
given in Table 6. 


Table 6, Geoid Comparison from Various Sources. 


Author 

Area 

Mean Diff. 
Auth. - Eapp 

Standard 

Deviation 

Marsh et al. 

West. Atlantic 

-3. 3 m 


Marsh et al. 

NE Pacific 

-3.2 


Marsh et al. 

SW of Aust. 

-2.1 


Brace (1977) 

5 areas 

4.6 
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The Brace undulations are based on Doppler satellite orbits and on a 
crossing arc adjustment process. 

The mean differences can be due to equatorial radius used as a reference 
(6378140 m for Marsh, and 6378135 m for Brace) and/or an uncorrected bias 
term. We will assume the bias has been correctly considered in Marsh while 
a correction of 2, 5 m is to be subtracted from the Brace data ( Anderle, 1979 ). 
The resultant mean differences then will imply an ideal equatorial radius as 
given in Table 7. 

Table 7. Equatorial Radius Implied by the Altimeter Geoid Comparisons. 


Comparison with: 

Equatorial Radius 

Marsh et al 
Brace 

6378137.1 m 
6378137.1 m 


The agreement of the two values is fortuitous. We must remember, however, 
that these values assume the orbit scale is correct. A formal accuracy assess- 
ment was not carried out but ±2 meters seems reasonable. 


The Computation of Mean Covariance Functions 

The estimation of the mean anomalies and mean undulations was to be 
done by the method of least squares collocation using the procedures described 
in Rapp (1977a). As part of these procedures it is necessary to determine point 
and mean covariance functions. As originally used in Rapp (ibid) the mean co- 
variance functions were determined from the numerical integration of point 
covariance fonctions. An alternative method, used earlier (Tscheming and Rapp, 
1974), by Schwarz (1976) and others, is to introduce the smoothing operator for 
degree H , jSjj (Meissl, 1971) into the series expressions for the needed covari- 
ance functions. The mean covariances needed were cov (A^,Ag), cov (N,Ag), 
and cov ( N, N) where an unbarred quantity represents a point value and a barred 
quantity represents a mean value. In our test computations we shall restrict the 
mean value to a l°x 1° anomaly block 

A point anomaly covariance function, given with respect to a reference 
field of potential coefficients to degree ^ ( ref) can be computed from (Tscheming 
and Rapp, 1974): ^ 

cov(Ag,Ag) = ^ c^ (cos ^) (1) • 

> . jJ=<^(ref)-!-i 

where Cjg are anomaly degree variances and s = (Ra/R)^ where Rb is the radius 
of the Bjerhammer sphere (internal to the earth) and R is the mean radius of the 
earth. One estimate for s is 0.999617 ( Tscheming and Rapp, 1974). To obtain 
the mean covariances we introduce : 
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( 2 ) 


cov (Ag,Ag) = y P^(cos 

^=^e(ref)+i 

The smoothing operator j3^ can be computed from the following (Rapp, 1977b) 

i8. = coti!^ ^(cos i/>o) ^ 3 ^ 

2 A(j&+1) 

where 0o is the radius of a cap that has the same area as the block (e. g. l°x 1°) 
being considered. As we are dealing with l°x 1° values the areas will be latitude 
dependent so that , in our application, will be latitude dependent. 

The other covariances needed can be found by the law of propagation of 
covariances (Moritz, 1972, 1978), In spherical approximation we have ; 

CO 

cov (N,Ag) = (R/G) ^ c^/(je-l) s^^ p^(cosi/)) (4) 

i=£(ref)+i 

CO 

cov (N,N) = (R/G)® ^ c_j/{^-l)®,%s^"‘^ p^(cosi^)) (5) 

Ji=^(ref) +1 

where G is an average value of gravity. Ih actual computations the summation 
to was replaced by a summation to degree 540 which is sufficiently high so 
that no significant change will take place by going to a higher degree. Tables of 
these mean covariances were constructed at 0?05 intervals. The original pro- 
gram was modified by removing the numerical integration procedure for the 
l°x 1° block computations and replacing it with table interpolation procedures. 

These methods were tested by carrying out predictions in two 5° equal 
area anomaly blocks both of which had a north latitude of -lO*^. The series co- 
variance functions were computed with a value for l°x 1° blocks at the equa- 
tor. Predictions were carried out using the numerical integration procedure and 
the directly tabulated mean covariance functions. The differences in the l°x 1 ° 
mean anomaly predictions was about ±0.6 mgals while the mean undulation 
differences were on the order of ±.lm. The estimated standard deviations 
were essentially unchanged. We thus conclude that this procedure could be 
applied successfully. 

The computer time savings were clearly seen in a run where the error 
covariance matrix (see later discussion) was also computed. The savings was 
14 secs ( on an IBM 370/168) or approximately 8% of the total time. This is a 
representative value only as the specific savii^s will depend on the number of 
given data points ( alternate undulations ). 

Although the savir^s..could be significant over repeated computations, the 
decision was made not to implement these procedures in our operational runs. 
One reason for this is that we felt that we would have to compute new covariances 
for different latitude areas because of the change in jS's with the area of the 
block. ( For example, for a l°x 1° at the equator 6 = 0° 564 while at latitude 
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60°, jS = 0?402. ) It could be that sufficient accuracy could be found by adopting 
some average values but this was not tried. Additional testing in this area 
is needed as well as the testing of series expression for the 5° equal area pre- 
dictions that were also carried out. 


The Error Covariance Matrix 


hi the computations performed previously (Rapp, 1977a) for mean anomalies 
and undulations, the l°x 1° values within a 5° equal area block were predicted 
from the same data set selected in and around the 5° block. This procedure had 
the advantage that only one matrix inversion was required; it did however have 
the disadvantage that the predictions of the l°x 1° blocks could be highly corre- 
lated. To specifically consider this question we considered the following ex- 
pression for the error-covariance matrix for the least squares collocation pro- 
cedure (Moritz, 1972); 


Ess ~ Css ~ Csx C CxB 


(6) 


where: 


Ess is the error covariance matrix of the signals being predicted; 

Css is the (physical) covariance matrix of the signals being predicted; 
Csx is the cross covariance matrix between the signals being predicted 
and the observations; 

C is the covariance matrix of the observations plus the noise covar- 
iance matrix; 


In our case the observations are the altimeter derived undulations while the 
signals being predicted are the l°x 1° mean gravity anomalies ( or undulations). 
The size of the Ess matrix will depend on the number of l°x 1° blocks within 
the 5° equal area block. Near the equator Ess would be a 25 x 25 matrix while 
at latitude 65° Ess would be a 55 x 55 matrix. The diagonal elements of Egs 
would be the square of the predicted accuracy of the individual blocks. Ess can 
also be converted into a correlation coefficient matrix. 


Equation ( 6) was evaluated for several cases of possible interest. In the 
first case only one known undulation value was used to estimate all the l°x 1° 
anomalies in a 5° block. The maximum correlation found was 0.53. When 88 
data points were used in this same 5° block, the maximum correlation coefficient 
was now . 13. 

Two additional 5° blocks were also investigated. One 5° block (at^n = 
-12? 5) had 391 data points used for the prediction. This would be considered a 
block with fairly dense altimeter data. In this case the largest correlation coef- 
ficient was 0.20. A second block, in the same latitude, having a less dense data 
set (261 points), was also considered, hi this case the largest correlation coef- 
ficient reached 0. 34. 
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We thus see that with increasing data sets ihe correlation between the 
blocks decreases. For most computations we will be doing, we would expect 
the correlation coefficient to be 0.3 or less. Because of this small magnitude 
we can regard most anomalies as almost independently determined. We there- 
fore choose not to compute the Ega matrix for our operational worik. 


Reference Model Error 


The prediction of mean gravity anomalies and geoid undulations has been 
carried out with respect to a reference gravitational field model defined by the 
GEM 9 potential coefficients (Letch et al. , 1977) taken to degree 20. Specifically 
the anomaly and its estimated accuracy is computed from the following equations 
(Rapp, 1977a, 1978a); 


where Ag 
h 


~ CghR ( CnhR + P) + (^) 

CggR - CgHR (ChhR + D)"^ ChgR (8) 

is the predicted free-air gravity anomaly with respect to 
ellipsoidal gravity field; 

is a column vector of the altimeter implied geoid undulations ; 


CghR is the row vector containing the covariance (referred to the 
reference field ) between the anomaly being predicted and 
the given geoid undulation; 

Chhfi is the square, symmetric matrix containing the covariances 
(referred to the reference field) between the given geoid 
undulations. If there are n h values being used this matrix 
is n X n; 

D is the error-covariance matrix of the given geoid -undulations 
which was taken to be a diagonal matrix whose elements 
corresponded to the square of the standard deviation of the 
altimeter measurement; 

Cggfi the expected mean square value (referred to the reference field) 
in a global sample, of the anomaly being predicted; 


mg the predicted standard deviation of the predicted anomaly; 

AgR , the gravity anomaly and geoid undulation implied by the reference 
Nr , hp set of potential coefficients. 


Similar equations can be written for the estimation of the mean undulations. 
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Ill our previous computations no consideration was given to the error in 
the final result caused by errors in the GEM 9 coefficients. One way suggested 
by Colombo (private communication, 1979) is to add to the covariance matrices, 
given with respect to the reference field, a component implied by the errors in 
the potential coefficients. This component can be obtained from equation (2), 

(4), and (5) by setting = 1 and replacing the Cjj values by the errors in the 
anomaly degree variances implied by the errors in the potential coefficients. 

Such a computation was done and prediction results for the i°x 1® anomalies 
and undulations in a 5® block were compared. We found changes in the predicted 
l°x 1° anomalies of only about ±0.2 mgals with the standard deviations increasing 
about 0.3 mgal (from about ±7 mgal). The undulation values and their accuracy 
changed no greater than 0.2 meter. 

We thus conclude this the errors in the GEM 9 coefficients, when used as 
a reference field, do not significantly contribute, to our final error estimate, 
provided the above procedure is correct. 


Anomaly Differences Implied by the New and Old Adjustments 

In our previous report we had computed 9995 l°x 1® anomalies. By Feb- 
ruary 1978 this number had increased to 12144 values, all values being based on 
the adjusted undulations from our first adjustment process. After completing the 
second adjustment of the altimeter data we computed a large number of l®x 1® 
anomalies to compare with our earlier estimates. Some statistics on these com- 
parisons are shown in Table 8. 


Table 8. Comparison of l®x 1° Mean Anomalies Computed from New 
and Old Adjustments of the Altimeter Data. 


Area* 

Number 

Mean Diff . 

RMS Diff. 
( mgals ) 

Max. Diff. 
(mgals) 

Calibration 

489 

0.2 

±4.9 

34.3 

South America 

1675 

0.2 

±2.3 

17.2 


* see Rapp (1977a) 


The root mean square differences are all smaller than the predicted accuracy of 
the l®x 1® anomalies and thus such differences are not of great concern. There 
is some concern over a few large discrepancies ( such as the 34 mgal difference). 
The reason for fliis is not clear although it may be related to a singly bad data 
point. 


The anomalies from the..new and old adjustment were also compared to ter- 
restrial anomalies. These comparisons are shown in Table 9, 
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Table 9. Comparison of l°x 1° Mean Anomalies from the New and Old 
Adjustment to Terrestrial Data. 


Area 

Old Adjustment 

New Adjustment 

Number 

RMS Diff. 

Number 

RMS Diff. 

Calibration 

260 

± 8. 9 mgals 

265 

±10. 7 mgals 

South America 

382 

16.5 

491 

17.7 

Alaska 

492 

9.5 

524 

10.5 . 

Phillippines 

527 

12.5 

194 

12.4 


It seems clear that the old adjustment yields slightly better anomaly values. 
Again considering the predicted anomaly accuracies these differences are not 
significant. However they are sufficient to have us accept in our final results, 
anomalies predicted from the old adjustment where such data exists in sufficient 
quantity to obtain reliable predictions. 


Effect of the Mass of the Atmosphere 

The computations for the gravity anomalies made through equation (7) have 
made no assumption on the attraction of the mass of the atmosphere. For points 
internal to a spherical shell comprising the atmosphere the attraction of the at- 
mosphere is zero. However in gravity anomaly computations the current pro- 
cedure is to include the mass of the atmosphere within the mass of the reference 
ellipsoid. Thus for comparisons of altimeter derived anomalies with terrestrial 
anomalies derived using a gravity formula based on a reference ellipsoid that 
contains the mass of the atmosphere, a small correction to the altimeter anomaly 
is needed. This correction can be found by first defining a gravity anomaly with 
respect to an ellipsoid for which the mass of the atmosphere is not included; 

hgE = g - 7 e ( 9 ) 

where E indicates the earth mass, g is observed gravity (properly reduced) 
and is normal gravity. The corresponding anomaly when the mass of the 
atmosphere is included in the reference ellipsoid is; 


AgE+A 

= g - yE+ A 

(10) 

Age+ A 

= Ag£ - (7E+A - Ve) 

(11) 


For points (or blocks ) located at a zero elevation (7 e+ a - ) is 0. 87 mgals. 

( lAG, 1971). Thus, before comparing the. altimeter derived anomalies to terres- 
trial anomalies that woiild be referred to the gravity formula of the Geodetic 
Reference System 1967, (for example) 0. 87 mgals should be subtracted from the 
altimeter derived anomaly. 
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Effect of Systematic Undulation Error on the Anomaly Predictions 

The ”geoid undulations” used in our estimation process may have systematic 
errors caused by the neglect of sea surface topography with respect to the geoid. 
To see the effect we carried out a computation of the l°x 1° anomalies and undu- 
lation in a 5° equal area block using altimeter derived undulations from which a 
constant 1 meter had been removed. Comparison with the same predictions 
without the systematic change revealed a systematic difference in the anomaly 
predictions of 2.5 mgals. This change is below the average accuracy of the 
predicted anomalies by a factor of 3. 

Subsequent comparisons of the l°x 1® predicted anomalies to terrestrial 
anomalies indicates a mean difference of 0.5 mgals. This difference could be 
due to sea surface topography effects that do not average to zero, to an error in 
the equatorial gravity of the Geodetic Reference System 1967, or it may be sta- 
tistically insignificant. 


Effect of Covariance Functions Used on the Predicted Anomalies and Undulations 


The predictions defined by equation (7) require certain covariance functions. 
In Rapp (1977a) the covariance functions were obtained from subroutine COVA 
(Tscheming and Rapp, 1974) with respect to a degree 20 field based on a certain 
anomaly degree variance model. The covariances used are called global covari- 
ances in the sense that they are based on parameters representative of the global 
gravity field. One argument is that it would be more reasonable to use covariance 
functions specifically designed for a given area. It thus seems appropriate to 
consider the sensitivity of our predictions to the covariances used in the prediction 
process. 

To do this we choose to work with three different covariance functions 
carrying out predictions in l°x 1° blocks within several 5° equal area blocks. 

The first covariance function group was that implied by the anomaly degree 
variance model described in Tscheming and Rapp (1974). This model was used 
to compute the needed covariances with respect to a degree 20 reference field. 
These covariances were the ones used in all of our previous computations (Rapp, 
1977a). 

The second covariance function was based on an anomaly degree variance 
model developed by Jekeli (1978) and designated as the IL model. This model 
used the same form of the anomaly degree variance model as used by Tscheming 
and Rapp but computed the model parameters enforcing a low horizontal anomaly 
variance (Co) considerably less than the T/R model, and an anomaly correlation 
length 4 (i.e. C(^) = Co/2) that was about twice that of the T/R model. 

The third covariance function group considered was that described by Jordan 
(1972) for local use. Jordan discussed a third-order Markov undulation covariance 
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model, from whicli an anomaly covariance model, and an anomaly-undulation 
cross covariance model could be derived. Specifically we have: 


C 


C(N,N) = Co(N) (l+^ 

C (Ag,Ag) = Co(Ag) ~ 


( 12 ) 

(13) 


(Ag,N) 


= 2/Co^Co(Ag) 



In these expressions Co(N) is the variance of the undulations and Co(Ag) is the 
variance of the anomalies in the area under consideration with respect to a refer- 
ence field which is the GEM 9 field in our ease. I„ is the Bessel function of the 
first kind of order n, and K„ is the Bessel function of the second kind of order 
n. (Although these functions are not defined at r = 0, definition is possible at a 
value of r sufficiently small to approximate zero.) The value of r is the dis- 
tance between the points under consideration and D is known as the character- 
istic distance. The values of Co(N) and Co(Ag) are related as follows ; 


yCo(N) ^ yCo(Ag) 
/3 D v/T go 


(15) 


where go is an average value of gravity ( 979. 8 gals ). 


We need to develop a procedure to determine the two independent parameters 
of this model. To do this we use the fact that the mean anomaly (or mean undialation) 
variance can be computed knowing the point covariance function ( Heiskanen and 
Moritz, 1967, p. 276). Thus, if we subdivide a given area into n^ subdivisions 
we can write 

n n n n 

var(^) ^ J C((Ag,Ag) Xi,yj,Xi’,yj') (16) 

1 j i' j' 

where x ,y are coordinates within the block. A similar equation can be written 
for var(N). 

Now suppose we take an area for which we known var (Ag) and var (N) with 
respect to some reference field. Then we can use these values to- determine Co{N) 
and Co(Ag) by solving (16) numerically. Given these two values the characteris- 
tic distance D can be found from (15). The process is an iterative one since some 
starting D value is needed, in our tests, convergence was obtained in an average 
of 5 iterations by stopping when the change in D was less than 1 meter. 


In our application we chose to find the local parameters Co(Ag), and D that 
would be characteristic of the l®x 1® anomalies and undulations_in specific 5^ 
equal area blocks. The data used for the computation of var(Ag) and var(N) 
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was the anomaly and undulations predictions with respect to the GEM 9 field 
when the predictions were done with the covariances from COVA using the T/R 
anomaly degree variance model. Given these two values we then computed 
Co(Ag), D and Co(N) which were then used to generate the covariances from 
equatioi^ (12), (13), and (14). 

We selected four 5° equal area blocks to work with. These blocks are 
numbered (see Figure 4) as follows: 1058 ( a block stradling the 90° East 
Ridge ), 1060 (a block in the Indian Ocean with mild anomaly variation with 
medium altimeter coverage ), 602 ( a block that contains the Puerto Rican 
Trench and has dense altimeter coverage ), and block 465 ( off the east coast 
of the United States with dense altimeter coverage ). 

The values of /Co(N) = CTn » /Co(Ag) = o^g , D , and the anomaly correla- 
tion distance ( i.e. the distance at which C(r) = Co(Ag)/e"^ ( l/e = . 3678. . .) ). 
These values ( except for D ) are also given, for comparison purposes, for the 
COVA covariances with the T/R anomaly degree variance model, and the 1 L 
model of Jekeli. 


Table 10. Parameters of Covariance Functions. 


Block/ 

Function 

mgals 

CTn 

meters 

D 

meters 

Anomaly Correl. Dist. 
(km) 

1058 

16.09 

1.28 

64042 

88 

1060 

11. 34 

3.00 

212365 

289 

602 

121. 88 

7.11 

46841 

64 

465 

26.55 

2.10 

63637 

87 

COVA 

38.98 

3.64 

- 

199 

IL 

31.97 

3.62 

- 

154 


The large Cg value for block 602 reflects the very large ( about -250 mgal ) 
anomalies associated with the Puerto Rican Trench. 

For blocks 1058, 1060, and 602 predictions were carried out with the 
three different covariance functions. For block 465 only the COVA with T/R 
model and the Jordan model were used. 

The predicted anomalies were then compared to corresponding terrestrial 
estimates to see if one covariance function gave better predictions than another. 
The root mean square ( RMS) differences, the RMS terrestrial anomaly standard 
deviation, and the RMS altimeter derived anomaly standard deviation are given 
in Table 11. 
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Table 11. Comparison of l°x 1° Anomalies Derived Using Different 
Covariance Functions to Terrestrial Estimates. 


Block Number 1058; R]VE Terrestrial Anomaly S.D. = ±10. Smgals 


Covariance 

RWB 

RMS Alt. Anomaly 


Difference (mgals) 

S.D. (mgals) 

COVA 

±12.7 

±7.6 

Jekeli IL 

12.7 

6.5 

Jordan 

12.0 

4.8 

Block Number 1060; RMS Terrestrial Anomaly S.D. = ±10. Smgals 

Covariance 

RMS 

RMS Alt. Anomaly 


Difference (mgals) 

S.D. (mgals) 

COVA 

±11.9 

±6.8 

Jekeli 1 L 

11.9 

5.6 

Jordan 

11.9 

2.4 

Block Number 602; RMS Terrestrial Anomaly S.D. =±11, Smgals 

Covariance 

RMS 

RMS Alt. Anomaly 


Difference (mgals) 

S.D. (mgals) 

COVA 

±12.2 

±6.3 

Jekeli 1 L 

12,4 

5.1 

Jordan 

12.1 

9.7 

Block Number 465; RMS Terrestrial Anomaly S.D. = ±13. Bmgals 

Covariance 

RMS 

RMS Alt. Anomaly 


Difference (mgals) 

S.D. (mgals) 

COVA 

±9.0 

±6.4 

Jordan 

9.6 

4.7 


From Table 11 we see that no significant improvement or difference in the pre- 
diction results can be seen from the use of different covariance fiinctions. The 
most significant change occurs in the predicted accuracy of the anomalies where 
the Jordan function displays accuracies more related to the variations of the 
anomaly field in a specific area. 

We have also compared the anomalies from the Jekeli 1 L function and the 
Jordan function to the anomalies from the COVA 1 L function. These comparisons 
are shown in Table 12. 
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Table 12. Comparison of liie i°x 1° Anomalies Derived from 
Two Alternate Covariance Functions to the values 
Obtained from the Covariances of COVA. 


Block 

Covariance 

used 

Mean Diff. 
(mgals) 

RMS Diff. 
(mgals) 

1058 

Jekeli IL 

0.0 

0.3 

1058 

Jordan 

0.1 

2.3 

1060 

Jekeli IL 

0.0 

0.1 

1060 

Jordan 

0.7 

3.0 

602 

Jekeli IL 

0.1 

0.4 

602 

Jordan 

-2.6 

7.9 

465 

Jordan 

-0.3 

1 

1.5 


We see that the comparison with the Jekeli IL fimction shows differences 
on the order of 0. 3 mgals while the differences with the Jordan function are 
higher, being about ±3 mgals in three blocks rising to 8 mgals in block 602. 

The large differences in this latter block occur for anomalies on the order of 
“280 mgals, the largest difference being 20 mgals. Five specific l°x 1° 
anomalies and undulations for blocks along the trench are shown in Table 13 
as computed using the three different covariance functions previously discussed. 


Table 13. I'^x 1° Anomalies and Undulations in the Area 
of the Puerto Rican Trench. 


Northwest 

Comer 


Undulations 

(meters) 


g 

Terr. 


IL 

Jordan 

COVA 

IL 

Jordan 

20 

292 

-244 ±24 

-231 ±6 

-231 ±5 

-243 ±11 

-62. 4 ±.2 

-62. 3±.2 

-62.7±.4 

20 

293 

-282 ±14 

-264 ±6 

-263 ±6 

-284 ±10 

-65. 6 ±.2 

-65. 6 ±.2 

“66.1 ±.4 

20 

294 

-205 ±22 

-223 ±7 

“222 ±6 

-239 ±12 

-63. 5±. 2 

-63.5±.2 

-64.0 ±.4 

20 

295 

-166 ±10 

-160 ±6 

-159 ±5 

-172 ±10 

-59. 7 ±.2 

-59.7±.2 

-60.0 ±.4 

20 

296 

-162 ± 6 

-148 ±7 


-149 ±6 

-161 ±13 

-59.1 ±.3 

-59.1±.2 

-59. 5 ±.5 


hi this section we have examined the variability of the prediction process as 
a function of the covariances used. There are some indications that the use of a 
tailored covariance function may give slightly better predictions than global func- 
tions but the evidence is marginal. On the other hand the standard deviations from 
the use of the Jordan function may be more realistic. In some cases these standard 
deviations are more and in some less than obtained from the global functions. In 
order to apply the local function we must have representative values for the anom- 
aly and undulation variances. In our application this was simple as predictions had 
already been carried out. If such predictions had not been done we could have used 
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the terrestrial data to obtain an anomaly variance and the results of some de- 
tailed geoid computations to obtain the undulation variance. If we were working 
in gravtmetrically unsurveyed areas ( such as the southern oceans ) we could 
select a global function. 

For our production estimation process we choose to work with the covari- 
ances implied by COVA with the T/R anomaly degree variance model. This was 
primarily done for consisting purposes as there did not seem to be sufficient 
evidence to suggest that significant gains would be obtained from using a tailored 
covariance function. Perhaps, with more accurate determination of oceanic 
mean anomalies from terrestrial data, a more definitive comparison and con- 
clusions could be made. 


Production Estimation of Mean Anomalies and Undulations 


The estimation of l®x 1° and 5° mean anomalies (and undulation) using 
equation (7) took place using the old adjustment where data was sufficiently dense, 
and using the new adjustment data in other areas. Many of the old adjustment 
anomalies were taken from Rapp (1977a). The covariance function used was 
obtained from COVA (Tscheming and Rapp, 1974) with respect to a degree 20 
field. Other specific details of the prediction process are described in Rapp 
(1977a). The values from the old and the new adjustment were merged together 
to form a combined data set. This data set contained 29479 l°x 1° blocks and 
957 5° equal area blocks, A number of the l°x 1° predictions were made in 
land areas and in ocean areas where the altimeter data was sparse. This was 
done only because of the manner chosen for the estimation of ^ the l°x 1° 
anomalies wi&in a 5° equal area block. The more reliable anomalies are those 
having an accuracy of ± 15 mgals or better. There are 27466 such values 
whose location is shown in Figure 11, (Asa matter of interest Figure 12 shows 
the 20599 values where the accuracy is ±8 mgals or better. ) A listing of the 5° 
equal area anomalies and undulations, referred to an ellipsoid whose flattening 
is 1/298.256, is given in the appendix. A tape containing the l°x 1° altimeter 
derived anomalies is also available. The predicted- anomalies have been com- 
pared to a terrestrial data set called "June 78 delete 424". This data set is that 
terrestrial field described in Rapp (1978b) less 424 anomalies that had a very 
large difference with the altimeter values. The net data set available for com- 
parison purposes contained 38981 values. This data set was also used to gener- 
ate a 5° equal area anomaly field that was used for the 5° block comparisons. 

The 5° comparisons' were made using only those 5° terrestrial anomalies where 
the terrestrial standard deviation was 10 mgals or less. The l°x 1° compari- 
sons were made only when the accuracy estimates for both anomalies were 15 
mgals or better. The resulte are given in Table 14. 
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is ± 15 mgals or Better. 











Table 14. Comparison of 5° Equal Area, and l°x 1°, Altimeter Derived 
and Terrestrial Mean Eree Air Anomalies. 


Description 

5° EA Value 

l^x 1® Value 

Mean Difference (GEM 9-Terr. ) 

0 . 3 mgals 

— 

Mean Difference (Alt. - Terr. ) 

0.7 


RMS Difference (GEM 9 - Terr. ) 

±8.9 

23.4 

RMS Difference (Alt. - Terr.) 

6.8 

11.8 

RMS Terrestrial Accuracy 

4.8 

10.9 

RMS Altimeter Accuracy 

2.7 

7.8 

RMS Terrestrial Anomaly 

15.2 

27.7 

Maximum Difference 

35.9 

63.8 

Number of Differences > [y [ mgals 

10 * 

7t 

Number of blocks compared 

767 

10086 


* y = 20, . t y = 40 


We see that the altimeter anomalies have a better agreement with the ter- 
restrial anomalies than the GEM 9 anomalies (computed from potential coeffi- 
cients to degree 20) as would be expected. The average accuracy of the 5° al- 
timeter anomalies is 3 mgals while it is 4 mgals for the terrestrial data used. 
The RMS difference between the 5° terrestrial and altimeter anomalies of 6. 8 
mgals is somewhat greater than would be expected if the terrestrial and altim- 
eter accuracies were correct. We finally note a very small (0.7 mgal) system- 
atic difference between the terrestrial and altimeter anomalies. 

The l°x 1° anomaly comparisons show a 11,8 mgals RMS differences 
between the altimeter and terrestrial data with only a 0.5 mgals systematic 
difference. This RMS difference is somewhat smaller than would be expected 
from the average accuracy estimates of the two data types which is about 11 
mgals for the terrestrial data and 8 mgals for the altimeter values. 

The l°x 1° mean geoid undulations estimated from the altimeter data were 
compared to the GEM 9 undulations (computed to degree 20) where we found a 
mean difference of 0.0 meters and a root mean square difference of ±2.7 m. 
This difference is consistent with the expected value of ±3.2 meters. 


Estimation in Small Areas 


The proceeding computations have described the estimation of mean anom- 
alies and mean undulations. It is of interest to examine the values to be obtained 
on the scale smaller than the 1° x 1° anomalies considered in this report. One 
attempt at this is described in Rapp (1978a, 1979) where point anomaly profiles 
were constructed across the Bonin Trench and around a sea mount in the Gulf 
of Alaska, hi the case of the trench we saw an anomaly change of 443 mgals in 
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118 km. This large change indicated that there was high frequency anomaly in- 
formation within the altimeter derived undulations. However we did not have 
at that time ship gravity data to compare with our altimeter 'derived anomalies. 

Recently Detrick and Watts (1979) have described some investigations near 
the Ninety-East Ridge, These studies pointed out a region of available ship gravi- 
ty data with bathymetric data that would provide a usefol test area. 

We first constructed an altimeter geoid in the area of the ridge which is 
shown as Figure 13. This map shows the undulation contours at 1 m intervals 
with the altimeter tracks plotted. In addition the 4000 m depth contour taken 
fromSclater and Fisher (1974) has been plotted to aide in the identification of the 
ridge crest. The depth of the actual crest varies approximately between 2000 m 
and 3500 meters. 

This map has been prepared using the altimeter data of the first adjustment 
using a prediction interval of 0® 5 x 0? 5. In addition a few bad data points (newly 
discovered) were removed from the data set. Consequently this geoid will differ 
somewhat from that global representation described earlier. 

For much of the geoid map there seems to be no specific association with 
the ridge although some dependence may be seen in the more southern parts. 

To specifically test the anomaly prediction process ship data was obtained 
from Watts (private communication). The tracks obtained corresponded to some 
of the profiles described in Detrick and Watts (1979). 

The point anomalies were predicted from the altimeter data at the same 
points as the existing ship gravity measurements. The prediction was carried 
out using one or two data selections and matrix inversions per profile. The 
results of these predictions are shown in Figures 14 and 15 where we have plotted 
the altimeter derived geoid, the altimeter derived anomalies, the ship detesmiined 
anomalies and the measured bathymetry. 

The profile shown in Figure 14 corresponds to the data used to obtain 
profile 90 E-9 in Detrick and Watts (1979). The starting (S) and ending (E) 
coordinates are; (ps = -3°19, Xg = 85?115, = -4?00, Xe = 93? 515. The pro- 

file shown in Figure 15 corresponds to the data used to obtain profile 90 E-12 in 
Detrick and Watts (1979). The starting and ending coordinates are; <pg = -16? 755 
Xs = 8? 3225, <Pe = -17?507, Xe - 91? 874. 

We can see that the altimeter anomalies follow quite closely the ship data 
except for the high frequency components. There is a clear correlation of the 
altimeter derived anomalies with the bathymetry. On the other hand the geoid 
seems to vary quite smoothly across the ridge, showing a small bump at the 
ridge crest. 


- 31 - 





RNOMflLT (MGRLS) 

- 51.00 - 28.00 - 5.00 18.00 41.00 64.00 


FIGURE I UNDULATION, ANOMALY, AND BATHYMETRY 
ACROSS 90 EAST RIDGE AT 4 > — 3.6° 


: O 
■ O 



0.00 25.00 50.00 75.00 100.00 125,00 15^.00 

DISTRNCE (KM) 


- 33 - 


RNQMRLT (MGRLS) 

- 51.00 - 28.00 - 5.00 18.00 41.00 64.00 


FIGURE IS. UNDULATION, ANOMALY, AND BATHYMETRY 
ACROSS 90 EAST RIDGE AT 4 > — 17° 




The R]VIS differences between the ship data and the altimeter data for the 
profile shown in Figures 14 and 15 was 12 mgals and 20 mgals respectively. The 
average predicted accuracy was 26 mgals. l£ the ship data is filtered to remove 
the high frequency components this agreement between the ship and the altimeter 
derived anomalies improves. For example, Eren (1979, private communication) 
has shown that if the components below a wavelength of 100 km are filtered out, 
the RMS differences for both profiles shown in Figure 14 and 15 become ± 11 
mgals. 

There seems to be two conclusions to be reached from these specific studies. 
First the altimeter derived anomalies agree quite well with the ship data 'consider- 
ing the accuracy of both data types. The altimeter anomalies seem to reflect the 
general bathymetry but they do not reflect the high frequency information seen by 
the ship measurements. More detailed analysis is needed to assess the accuracy 
of both data types by wavelength. In addition we need to look at the use of more 
detailed altimeter data as opposed to using from averages which represent an 
immediate averagir^ and loss of high frequency information. 

Second we see that the geoid undulations only slightly reflect the variations 
going over the ridge. The significant variation is seen much more in examining 
the anomaly data. 


Summary and Conclusions 


This report has described the processing and analysis of Geos-3 data re- 
leased for general use by March 1978. The data was edited and adjusted to re- 
move bad data and to remove orbit and altimeter bias terms. The adjustment 
was carried out first in a primary set involving 700 arcs chosen for their global 
distribution, and 4 regional adjustments. The crossover discrepancies after the 
adjustment averaged about ±55 cm. 

The resultant data waS' used to prepare a global sea surface map ( or approx- 
imately the geoid) with accuracy estimates. Comparisons of these maps with 
similar data produced using precise orbits indicated random differences on the 
order of 1 meter with a systematic difference implying an equatorial radius of 
6378137 meters. Computations were also done to investigate the accuracy of 
the altimeter geoid by wavelength. Comparisons with the undulation implied by 
the GEM 9 potential coefficients showed differences of about 75 cm at 5 wave- 
lengths from 13010 km to 2602 km. 

Several investigations were also carried out to improve our mean anomaly 
and mean undulation computation procedures. Specifically we considered the use 
of the smoothing operator in the series covariance expressions instead of numer- 
ical integration procedures for mean covariance computations. Although the pro- 
cedure could save some computer time without any significant loss of accuracy, 
the original numerical integration procedure was retained for logistical reasons. 
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We also considered the correlation between various l°x predictions within a 
5° equal area block. We found that the average correlation coefficient between 
adjacent blocks was about 0.2 so that the blocks could be considered to have 
been determined independently. 

The altimeter data was then used to extend our 5 ° and l°x 1° mean anom- 
aly and undulation computation. For some unknown reason the data from our 
original adjustment in 1977 gave somewhat better anomalies than the new ad- 
justment. Consequently in those geographic areas where the old adjustment 
data was adequate we used that data for anomaly estimation. In the other areas 
the new adjustment data was used. A total of 29479 l°x 1® blocks and 957 5° 
equal area blocks were estimated. Of these l°x 1° values 27466 had predicted 
standard deviations of ± 15 mgals or better. Comparisons were made with the 
terrestrial anomaly data where we found differences on the order of ±12 mgals 
for the l°x 1*^ data and ±7 mgals for the 5° data, A good part of this difference 
is due to the errors in the terrestrial anomaly field. 

The anomaly and undulation predictions were done with a global covariance 
function referred to a degree 20 reference field. Additional tests were carried 
out with a different global covariance function and a local covariance function 
derived for special areas. No significant changes were seen in the predictions 
from the two global field models. Root mean square differences on the order of 
3 mgals were found between the global and local covariance fimction but no sig- 
nificant improvement in the prediction process was seen based on comparisons 
with actual data. We did see chaises in the predicted standard deviations up to 
65% when using the local model instead of the global covariance model. The 
difficulty in applying a local model lies in a need for an adequate knowl edge of 
the residual anomaly and undulation field in the area. 

Point anomaly and undulation values were also computed along ship tracks 
that crossed the Ninety-East Ridge. The resultant anomaly values were compared 
to ship data where root mean square differences were 12 and 20 mgals for the two 
tracks considered. Plots of this data with the bathymetric indicate the correla- 
tion of the altimeter anomalies with the bathymetry on a regional basis as opposed 
to a more local correlation seen in the ship gravity data. Part of this difference 
may be due to our use of altimeter averages over 14 to 20 km swaths. However, 
we clearly can see the potential for improved knowledge of the anomaly field in 
small areas using the altimeter data. 

What additional things need to be done ? First the data needs to be reexamined 
to try to delete some additional bad data that has shown up. We should try to get 
additional data to fill in those areas for which data is non-existent or sparse. Then 
a readjustment could be done and the anomalies and undulations computed using 
localized covariance functions. One should also examine the use of non-frame 
averages going back, perhaps to the original data, or 1 second averages. The 
information, by wavelength, in the altimeter data should be studied to examine the 
accuracy of our anomaly and undulation results. And, of course, attempts have 
been made, and could continue to define sea surface topography effects. 
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Finally, it sould be clear that this Geos-3 altimeter data has enabled a 
significant improvement in our knowledge of the earth’s gravity field. It's use 
in local areas and in global computations has added significantly to geodetic and 
geophysical research. This data and Seasat-1 data will provide a data source 
for additional research into the gravity field at sea and indirectly on a global 
basis. 
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Appendix 


This appendix contains a listing of the 5° equal area free-air anomalies 
and undulations as derived from Geos-3 altimeter data. These values refer to 
an ellipsoid whose flattening is 1/298.257 and whose equatorial radius is theo- 
retically xmknown, being actually the equatorial radius of the general terrestrial 
ellipsoid. 
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0,2 

842 

0 

-5 

75 

70 

-39 . 9 . 

2.6 

- 89.6 

0 . 1 

843 

0 

-5 

80 

75 

- 46 . 1 

2.5 

- 96.6 

0.1 

844 

0 

-3 

85 

80 

- 42.3 

2.0 

- 91.8 

0 , 1 

845 

© 

-5 

9 © 

85 

- 23.9 

2.5 

- 74.0 

0 . 1 

846 

0 

-3 

95 

90 

- 8.5 

2.6 

- 50 , 1 

0 , 1 

847 

0 

-5 

100 

95 

10.5 

2.6 

- 24.0 

0 . 1 

848 

0 

-5 

105 

100 

35.2 

4,0 

2. 1 

0.6 

849 

0 

-5 

110 

105 

38.0 

2*8 

23.4 

0 . 1 

853 

0 

~5 

130 

125 

15,5 

2.7 

60.4 

0 . 1 

854 

0 

-5 

135 

130 

81.2 

2.6 

67.6 

0.2 

855 

0 

-5 

140 

135 

29.0 

3.4 

73 . 1 

0.5 

856 

0 

-5 

145 

140 

87.2 

3. 1 

76.3 

0*4 

857 

0 

-5 

160 

145 

44,7 

2.0 

75 . 5 

0 . 1 

858 

0 

-5 

155 

15 © 

21.0 

2.5 

65.9 

0 . 1 

859 

0 

-5 

160 

155 

20 . 1 

2.0 

59.6 

0 , 1 

860 

0 

-5 

165 

160 

3 . 1 

2.6 

48. 3 

0 . 1 

863 

0 

-5 

180 

175 

- 6.7 

2.0 

25.5 

0 . 1 

S 64 

0 

-5 

185 

18 © 

- 8,2 

2.4 

20.4 

0 . 1 

865 

© 

-5 

190 

185 

- 3.9 

2.4 

17.6 

0 , 1 

866 

0 

-5 

195 

190 

- 2.4 

e^txj 

15.9 

0 , 1 

867 

0 

-5 

200 

195 

3.8 

2.4 

15.9 

0 . 1 

868 

0 

-5 

205 

200 

10.7 

lil ♦ ^ 

16 . © 

0. 1 

869 

0 

-5 

210 

205 

14.7 

2.5 

14.3 

0 . 1 

870 

0 

-5 

215 

210 

17.0 

O r? 

t. 

11.1 

0 . i 

877 

© 

-5 

250 

245 

- 2.6 

2.5 

- 17.6 

0 . 1 

878 

© 

-5 

255 

25 © 


i:il m 

-IS . 7 

0. 1 

879 

0 

-5 

260 

255 

- 7.5 

2.6 

- 14.8 

0 . 1 

880 

0 

-S 

265 

20 © 

- 6.7 


- 12.7 

0 . i 

881 

0 

-5 

270 

265 

-S .0 

2. 5 

- 9.4 

0 . 1 

882 

0 

— o 

275 

270 

— 4.4 


- 2.9 

0. 1 

883 

0 

-3 

280 

275 

2.2 

2.9 

6 . 0 

0.2 

893 

© 

-5 

33© 

325 

- 3,7 

2.5 

- 4 . 1 

0 . 1 

894 

0 

-5 

835 

330 

0.4 

2,5 

2.2 

0 . 1 

895 

0 

-5 

84 © 

336 

- 2.5 


7,4 

0 . 1 

396 

0 

-5 

345 

340 

5,8 

2.0 

14.7 

0 . 1 

897 

0 

-5 

360 

S 45 

11,7 


20.3 

0 , 1 



SEQ 

LAT 

LOKG 

AlfOII 

S .». 

DKD 

S . D . 

SEa 

LAT 

LONG 

ANOM 

S . D . 

UHD 

S . D . 

898 

0 

-5 

355 

350 

4.3 

2.5 

19.7 

0. 1 

899 

0 

-5 

0 

355 

4.1 

2.5 

18.7 

0 . 1 

900 

-5 

-10 

5 

0 

- 5.2 

2.6 

11.9 

0 . 1 

901 

-5 

- 1 © 

10 

5 

- 7 . 1 

2.8 

9.6 

0.2 

907 

-5 

-10 

41 

35 

“ 9.9 

5.8 

- 27.5 

1.3 

908 

-5 

-10 

46 

41 

- 26.1 

2.4 

- 35,1 

0 . 1 

909 

-5 

-10 

51 

46 

- 18.6 

2.5 

- 87,7 

0 . 1 

910 

-5 

-10 

56 

51 

- 4.3 

Sis 

- 07.7 

0 , 1 

911 

-•5 

-10 

61 

56 

0.0 

2.7 

- 42.8 

0 . 1 

912 

-5 

-10 

66 

61 

- 13.4 

2.3 

- 54.6 

0.2 

913 

-5 

-10 

71 

66 

- 21,1 

2.5 

- 65 . 9 

0 . 1 

914 

-5 

- 1 © 

76 

71 

- 26.8 

2.7 

- 77 . 0 

0. 1 

915 

-5 

-10 

81 

76 

- 39 . 1 

2,8 

- 83.4 

0,2 

916 

-S 

- 1 © 

86 

81 

- 32.3 

2.5 

- 79.0 

0 , 1 

917 

-5 

-10 

91 

86 

- 24.9 

2.5 

- 67.0 

0. 1 

918 

-5 

-10 

96 

91 

- 16. 1 

2.6 

- 49.7 

0 . 1 

919 

-5 

-10 

101 

96 

8.4 

2.3 

- 27.0 

0. 1 

920 

-5 

-10 

106 

101 

1.7 

2.7 

- 7.7 

0.2 

921 

-5 

-10 

112 

106 

31.8 

2.6 

18.5 

0 . 1 

922 

-5 

-10 

117 

112 

33.8 

2.5 

82 . 0 

0 . 1 

923 

-5 

- 1 © 

122 

117 

38.7 

2.9 

41.7 

0.2 

925 

-5 

- 1 © 

132 

127 

- 23,6 

2.5 

47.5 

0.1 

926 

-S 

-10 

137 

132 

21.5 

2.3 

60.0 

0 . 1 

929 

-5 

-10 

152 

147 ■ 

29 . 1 

2.7 

72.8 

©. 1 

930 

-5 

-10 

157 

152 

20 . 1 

2.5 

67.6 

0 . 1 

981 

-5 

- 1 © 

162 

157 

16.5 

2.6 

62.1 

0 , 3 

934 

-5 

-10 

177 

172 

- 5.2 

2.4 

87 . 1 

0 . 1 

985 

-5 

-10 

133 

177 

- 10.9 

2.5 

35.3 

0 . 1 

936 

-S 

-10 

183 

183 

- 11.7 

2.4 

22.2 

0 . 1 

937 

“5 

- 1 © 

193 

133 

- 8.0 

2.4 

13.3 

0 . 1 

938 

-5 

-10 

193 

193 

1.3 

2.5 

16 . 4 

0 , 1 

939 

“5 

-10 

203 

19 S 

1.4 

0.4 

13.5 

0.1 

940 

-5 

-10 

208 

203 

5.3 

2.5 

11.4 

0 , 1 

941 

-5 

-10 

21 S 

20 S 

6.2 

2,5 

8.4 

0 . 1 

946 

-5 

-10 

233 

233 

0.0 

2.8 

- 11.7 

0 . 1 

947 

-5 

-10 

243 

23 S 

- 1.6 

2.5 

- 1 S .0 

0 . 1 

949 

-5 

-10 

254 

248 

4.6 

2.6 

- 12.4 

0 . 1 

950 

-5 

- 1 © 

259 

254 

- 3. 1 

2.6 

- 11.6 

0 . 1 

951 

-5 

-)0 

264 

259 

- 8.8 

2.5 

- 12.7 

0.1 

952 

-5 

-10 

269 

264 

- 10. 1 

2.5 

- ii .3 

0 . 1 

953 

-5 

-I© 

274 

269 

- 11 . 1 

2.6 

- 7 . 1 

0.1 

954 

-5 

-10 

279 

274 

- U .4 

2.6 

- 0.2 

©. 1 

964 

-s 

-10 

330 

325 

- 10.4 

2.7 

- 7.3 

0 . 1 

966 

-5 

- 1 © 

335 

SS 0 

- 10.6 

2.7 

- 4.4 

0.2 

966 

-5 

-10 

340 

335 

- 8.4 

2.7 

1.3 

0 . 1 

967 

-5 

-10 

345 

340 

4.4 

2.5 

10.4 

0. 1 

968 

-5 

-10 

350 

345 

9 . 1 

2.5 

16.0 

0 . 1 

969 

-5 

-IQ 

356 

350 

Q .6 

2.4 

16.2 

1 

970 

“5 

-10 

0 

355 

6.6 

2.5 

16.0 

0 . 1 

971 

- 1 © 

-15 

5 

0 

- 3.2 

2.5 

12.2 

0. 1 

972 

-10 

-15 

10 

5 

0 . 1 

2 .S 

14 . 1 

0.2 

973 

- 1 © 

-15 

35 

10 

13.9 

3,6 

17.8 

0.6 

979 

- 1 © 

-15 

4.6 

41 

- 14 . 1 

2.5 

- 25,4 

0 . 1 

980 

- 1 © 

-15 

5 1 

46 

2.0 

2.6 

- 23.8 

0. 1 

981 

-10 

-15 

57 

51 

- 1 . 1 

2.6 

- 32.5 

0. 1 

932 

-10 

-15 

62 

57 

3.3 

2.6 

- 31.6 

0 , 1 

983 

-10 

-15 

67 

62 

- 1.8 

2.6 

- 39.2 

0,2 

984 

-10 

-15 

72 

67 

- 4 .B 

2,5 

- 49.0 

0. 1 

985 

- 1 © 

-15 

77 

72 

- 21.3 

2.5 

- 60.7 

0 . 1 

986 

-10 

-ts 

32 

77 

- 21 .© 

2.5 

— 65 . 4 

0.1 

987 

- 1 © 

-15 

87 

82 

- 24.3 

2.5 

- 65.7 

0 , 1 

988 

-10 

-15 

93 

87 

- 00.6 

2.6 

- 71.0 

0. 1 

989 

-10 

-15 

98 

93 

- 13.7 

2. 5 

— 45 . 3 

0 , 1 

990 

-10 

-15 

103 

93 

- 9.2 

2,5 

- 30,3 

0. 1 

991 

-10 

-15 

108 

103 

1.2 

2.5 

- 14.5 

0.1 

992 

-10 

-15 

113 

lOS 

- 1 .© 

2.5 

“ 1.6 

0,1 

993 

-10 

-15 

118 

110 

“ 14.9 

2,4 

9.2 

0 . 1 

994 

-10 

-IS 

123 

IIS 

- 2.2 

2.6 

24.2 

0 . 1 

995 

-10 

-15 

129 

123 

27.0 

2.8 

47.7 

0 . 1 

997 

-10 

-15 

139 

134 

23. 7 

2.9 

60, 4 

0 .S 

998 

- 1 © 

-15 

144 

139 

21,1 

0 . 1 

65 . 4 

0.8 

100 © 

-10 

-15 

1 S 4 

149 

U.5 

2.5 

67 . S 

0 . 1 

1003 

-10 

-15 

170 

165 

11,8 

2.5 

53.4 

0. 1 

1004 

-10 

-IS 

175 

170 

22.3 

2.5 

ns . 6 

0.1 

1005 

-10 

-15 

180 

175 

12.2 

2.5 

46.8 

0 . 1 

1006 

-10 

-15 

135 

18 © 

10.5 

2.S 

33,5 

0 , 1 

1007 

-10 

-15 

190 

185 

3.6 

2.3 

28.6 

0 . 1 

1003 

-10 

-15 

195 

19 © 

- 2.2 

2.5 

19.5 

0 , 1 

1009 

-10 

-15 

201 

195 

0.3 

2.5 

16.3 

0.1 

1010 

- 1 © 

-15 

206 

201 

- 1.0 

2.5 

8.9 

0.1 

1011 

-10 

-15 

211 

206 

- 1.0 

2.6 

5 . 4 

0 . 1 

1012 

-10 

-15 

216 

211 

“ 3.7 

2.6 

0.4 

0,1 

1015 

-10 

-15 

231 

226 

1.4 

2.6 

-3.3 

0.1 

1016 

- 1 © 

-IS 

£ 

281 

2.9 

2,6 

- 10,9 

0 . 1 

1017 

-10 

-15 

242 

237 

6 i 1 

2.5 

-7,7 

0.1 

1018 

- 1 © 

-15 

247 

242 

10 . S 

2.6 

- 5.6 

0.1 

1019 

-10 

-15 

252 

247 

11.0 

2 .S 

- 4.7 

0 . 1 

1020 

-10 

-15 

257 

252 

6,9 

2.I-; 

- 6.3 

0.1 

1021 

-10 

-15 

262 

257 

- 3.0 

2.5 

-O.S 

0 . 1 

1022 

-10 

-15 

267 

262 

- 1.4 

2 . 4 . 

-a. 3 

0 . i 

1023 

-10 

-15 

273 

267 

“ 6.2 

2.6 

-3.0 

0 . 1 

1024 

-10 

-15 

278 

273 

- 10.0 

2.5 

- 1.6 

0 . 1 

1025 

-10 

-15 

283 

278 

- 16,7 

2,6 

3.2 

0 . 1 

1034 

- 1 © 

-15 

329 

324 

- 7,4 

2.5 

- 7.9 

0, 1 

1035 

- 1 © 

-15 

334 

329 

- 10,5 

2.5 

- 6,4 

0 . 1 

1036 

- 1 © 

-15 

339 

SS 4 

- 10.1 

2.6 

- 2.7 

0 . 1 

1037 

-10 

-15 

845 

3 S 9 

-3.0 

2.8 

6.3 

0 . 1 

1038 

-10 

-15 

350 

S 45 

1.5 

2.5 

H .2 

0 . t 

1039 

-10 

-15 

355 

350 

1 , 1 

2.4 

13.0 

0 . 1 

1040 

- 1 © 

-15 

0 

355 

— 1.2 

2.6 

12.9 

0,1 

104 ) 

-15 

-20 

5 

0 

- 0.8 

2.5 

io . 3 

0.1 

1042 

-16 

-20 

(© 

5 


2-4 

17.7 

0. 1 

1048 

“15 

-20 

42 

37 

- 13.7 

3.0 

- 12.6 

0.3 

1049 

-IS 

-20 

47 

42 

0.9 

3.4 

- 12.3 

0.6 

1050 

-15 

-20 

52 

47 

14.9 

3.5 

- 10.5 

0.5 

1051 

-15 

-20 

57 

52 

6,9 

2.6 

- 14.3 

0.2 

1052 

-15 

-20 

63 

57 

14.0 

2.7 

- 20.5 

0. 1 

1053 

-16 

- 2 © 

60 

63 

13 , 1 

2.7 

- 21.7 

0.2 



sEa 

LAT 

LONG 

AROM 

S . D . 

UNO 

S . D . 

SEa 

LAT 

LONG 

ANOM 

S . D . 

UND 

S . D . 

1054 

-15 

-20 

73 

68 

4.6 

2.5 

- 31.1 

0. 1 

1055 

-15 

-20 

73 

73 

- 4.7 

2.5 

- 41.3 

0. 1 

1056 

-15 

-20 

83 

78 

- 11.2 

2.6 

- 48.6 

0. 1 

1057 

-15 

- 2 © 

89 

88 

- 18.0 

2,7 

- 62.6 

0. 1 

105 & 

-15 

-20 

94 

89 

- 22.4 

2.5 

- 51.4 

0. 1 

1059 

-15 

-20 

99 

94 

- 22.6 

2.5 

- 46.4 

0. 1 

1060 

-15 

-20 

104 

99 

- 26.0 

2.4 

- 89.0 

0.1 

1061 

-15 

-20 

110 

104 

- 25.9 

2,5 

- S 3 ,3 

0. 1 

1062 

-16 

-20 

115 

110 

- 17.7 

2.4 

- 12.8 

0. 1 

1063 

-15 

-20 

120 

ns 

- 13.4 

2.7 

1.4 

0. 1 

1064 

-15 

-20 

125 

120 

13.4 

0.7 

19.0 

0.7 

1069 

-15 

-20 

151 

146 

IS . 9 

2.9 

61.4 

0.2 

1070 

-15 

-20 

157 

151 

19.5 

2.0 

70.7 

0. 1 

1072 

-15 

-20 

167 

162 

17.9 

2.5 

00.4 

0 . 1 

1073 

-15 

-20 

172 

167 

31.4 

2.5 

63.7 

0 . 1 

1074 

-15 

- 2 © 

177 

172 

82.4 

2.5 

61.2 

0 . 1 

1075 

-15 

-20 

183 

177 

25.9 

2.5 

62.0 

0 . 1 

1076 

-15 

-20 

183 

ISS 

22.5 

2.5 

41.6 

0. 1 

1077 

-15 

-20 

193 

188 

6.7 

2.4 

25.4 

0 . 1 

1078 

-15 

-20 

198 

19 S 

—3 • 6 

2.4 

15.4 

0 . 1 

1070 

-15 

-20 

203 

198 

0.5 

2.0 

10.2 

0 . 1 

1080 

-15 

-20 

209 

203 

5 . 1 

2.5 

7.0 

0 . 1 

1081 

-15 

-20 

214 

209 

4.0 

2.0 

1 . 1 

0 . 1 

1084 

-15 

- 2 © 

236 

224 

- 5 .S 

2.5 

- 12.3 

0 . 1 

1085 

-15 

-20 

235 

230 

0 .S 

2.4 

- 9.1 

0 . 1 

1086 

-15 

-20 

240 

2 S 5 

6.5 

2.5 

“ 6.4 

0 . i 

1087 

-15 

-20 

245 

240 

11.1 

2.0 

- 3.4 

0 . 1 

1088 

-16 

-20 

25 © 

245 

12.7 

2.5 

- 1.8 

0 . 1 

1089 

-15 

-20 

256 

250 

11.2 

2.7 

- 3.7 

0.2 

1090 

-15 

-20 

26 J 

256 

4.0 

2.5 

“ 4,7 

0 . 1 

1091 

-15 

- 2 © 

266 

261 

2.4 

2.5 

- 4.6 

0 . 1 

1092 

-16 

-20 

271 

266 

- 0.3 

2.7 

- 3.4 

0.2 

1093 

-15 

-20 

277 

271 

- 5.0 

2.6 

- i .0 

0 . 1 

1094 

-IS 

- 2 © 

282 

27 ? 

- 0.7 

2.5 

7.7 

0. 1 

1095 

-15 

-20 

287 

282 

1.0 

2.7 

18.0 

9.2 

1103 

-15 

-20 

329 

323 

- 2.9 

2.7 

- 7.0 

0.2 

1104 

-15 

-20 

334 

329 

- 7.4 

2.5 

- 5. 1 

0 . 1 

1105 

-15 

-20 

339 

SS 4 

- 13.3 

• f 

“ 3.8 

0. 1 

1106 

-15 

-20 

344 

339 

- 2.0 

2.7 

3.6 

0.2 

1107 

-15 

-20 

350 


3.5 

2.5 

12 . 1 

©. 1 

1108 

-)5 

-20 

355 

350 

1.1 

2.6 

12.8 

0 . 1 

1110 

- 2 © 

-25 

5 


5.9 

2.6 

21 . 1 

0 . 1 

nil 

-20 

-25 

11 

5 

12.4 

2.0 

29.7 

0 . 1 

1116 

-20 

-25 

33 

iS'^4 

- 5 .S 

4.0 

4.7 

0.7 

1117 

-20 

-25 

43 

38 

- 1.3 

2.. 5 

- 0.6 

0 . 1 

1119 

-20 

-25 

54 

48 

- 4 . 1 

2.7 

- 3.8 

0 . 1 

1120 

-20 

-25 

59 

54 

10 . 1 

2.6 

- 2.8 

0 . 1 

1121 

- 2 © 

-25 

64 

59 

12.5 

2.6 

“4 . 3 

0 . 1 

1122 

-20 

-25 

7 © 

64 

21.4 

2.6 

- 8.9 

0 . 1 

1120 

-20 

-25 

75 

70 

11.2 

2.5 

“ 16.5 

0 . 1 

1124 

-20 

-25 

81 

75 

0.0 

2.0 

- 33.5 

0 . 1 

1125 

-20 

-25 

86 

81 

- 6.3 

2.6 

- 36.2 

0.2 

1126 

-20 

-25 

91 

86 

- 7. 1 

2.5 

- 40.2 

0 . 1 

1127 

-20 

-25 

97 

91 

- 20.3 

2,5 

- 63.0 

0 . 1 

1128 

-20 

-25 

102 

97 

- 25.8 

2.4 

- 43.3 

0 . 1 

1129 

-20 

-25 

107 

102 

- 16.5 

2.4 

- 35 . 9 

0. 1 

1130 

-20 

-25 

113 

107 

- 22.5 

2.5 

- 30 . 1 

0. 1 

1131 

-20 

-25 

118 
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